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ABSTRACT 

, Context. We present the results of a long M87 monitoring campaign in very high energy y-rays with the MAGIC-I Cherenkov telescope. 
Aims. We aim to model the persistent non-thermal jet emission by monitoring and characterizing the very high energy y-ray emission of M87 
\l ' during a low state. 

Methods. A total of 150 h of data were taken between 2005 and 2007 with the single MAGIC-I telescope, out of which 128.6 h survived the data 
quality selection. We also collected data in the X-ray and Ferrai-LAT bands from the literature (partially contemporaneous). 
■ Results. No flaring activity was found during the campaign. The source was found to be in a persistent low-emission state, which was at a 
' confidence level of la. We present the spectrum between 100 GeV and 2TeV, which is consistent with a simple power law with a photon index 
\ T = 2.21 ± 0.21 and a flux normalization at 300 GeV of (7.7 ± 1.3) x 10~ 8 TeV" 1 s" 1 nr 2 . The extrapolation of the MAGIC spectrum into the GeV 
I , energy range matches the previously published Fermi-LAT spectrum well, covering a combined energy range of four orders of magnitude with 
• • ■ the same spectral index. We model the broad band energy spectrum with a spine layer model, which can satisfactorily describe our data. 
> 

• i-H ■ Key words, gamma-rays: galaxies — individual: M87 

X: 

a- 

1 . Introduction SMBH are able to energise and collimate the flow of matter in 

the form of relativistic jets, which extend to distances of kpc 
Current astrophysical observations suggest that the centre of al- up t0 Mpc from the centre. The exact mechanisms of this pro- 
most every galaxy harbours a super-massive black hole (SMBH) cess> as well as the source of energy that accelerates particles to 
with a mass of the order of 10 6 - 1O 1() M . Depending on the re lativistic velocities, are not yet entirely understood (see for in- 
evolutionary stage and the type of the galaxy, the amount of in- stance the Blandford-Payne scenario Blandford & Payne 1982, 
terstellar matter in the vicinity of the central black hole differs or the Blandford-Znajek process Blandford & Znajek 1977). 
significantly. In about 1% of the galaxies the gravitational en- 
ergy released by matter falling onto the central SMBH is suf- Once ejected through the jet, the particles interact with each 
ficiently high to make the nucleus more luminous than the en- other and the surrounding magnetic and radiation fields produc- 
tive galaxy, hence the name of active galactic nuclei (AGN, see ing photons with energies up to TeV (Marashi et al. 1992; Bloom 
Urry & Padovani 1995; Urry 2003; Vellieux 2003). In a small & Marscher 1996). The collimation of the jet and the relativistic 
group of AGN, which are called radio-loud, mechanisms (me- boosting of the emitted radiation imply that relativistic jets ap- 
diated by the magnetic field) acting close to the vicinity of the pear brightest when they are inclined by a small angle towards 
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the observer, as is the case of the sources known as blazars. 
Owing to their large distances and the very low inclination of 
their jets to the line of sight, it is difficult to study the structure 
of these sources. However, observation of some nearby radio- 
loud AGN that have a higher jet inclination angle (such as M87) 
give us the opportunity to study jets in greater detail, since the 
angular resolutions of radio to X-ray observatories is sufficient 
to resolve their structures (Wilson & Yang 2002). 

M87 is a giant elliptical radio galaxy at the centre of the 
Virgo cluster at a distance of 16.7 Mpc (Mei et al. 2007). It is 
harbouring a SMBH with a mass of (6.4+0.5)x 1O 9 M (Gebhardt 
& Thomas 2009). The prominent jet of M87, first discovered in 
the visual observation of Curtis (1918), is inclined by 10° - 45° 
from our line of sight (Biretta, Sparks & Macchetto 1999; Ly, 
Walker & Junor 2007). The vicinity of M87 makes it possible to 
study different parts of the jet in detail. 

The emission of the compact nucleus as seen in GHz fre- 
quencies is variable and has a dimension of several hundred 
Schwarzschild radii. Throughout the inner section of the jet sev- 
eral brighter spots, so-called "knots", are visible, often charac- 
terized by superluminal motion with apparent speeds reaching 
~6 c (Biretta, Zhou, & Owen 1995; Biretta, Sparks & Macchetto 
1999). The one closest to the nucleus (60 pc), known as HST- 
1, contains blobs with superluminal motion, as revealed by ra- 
dio VLBI observations (Cheung et al. 2007), and shows high 
variability of its optical and X-ray emission (e.g. Perlman et al. 
2003, Harris et al. 2009). A closer inspection also indicates the 
existence of a counter-jet (Biretta, Sparks & Macchetto 1999). 
Recently, Had a et al. (20 1 1 )| used the core-shift effect detected 
in Very Long Baseline Array (VLBA) multi-frequency data be- 
tween 2 and 43 GHz to estimate that the central engine of M87 
lies within a distance of 14-23 Schwarzschild radii of the radio 
core at 43 GHz. 

The first hint of very high energy (VHE, E>100GeV) y- 
ray emission reported by the HEGRA collaboration (Aharonian 
et al. 2003) triggered extensive observations by the next gen- 
eration Imaging Atmospheric Cherenkov telescopes (IACTs): 
The H.E.S.S. collaboration firmly established M87 as an emitter 
above 730 GeV (Aharonian et al. 2006) and revealed flux vari- 
ability on time-scales of two days, suggesting that the emission 
region of the y-rays is very compact with a dimension similar 
to the Schwarzschild radius of the central black hole. VERITAS 
also detected VHE y-radiation from M87 in 2007 (Acciari et al. 
2008) above 250 GeV and subsequently monitored the source 
during the following years (see e.g. Acciari et al. 2010). The first 
reported detection of y-ray emission from M87 by the Major 
Atmospheric Gamma-Ray Imaging Cherenkov (MAGIC) col- 
laboration used data taken in 2008 during a flaring state (Albert 
et al. 2008b). These data were collected during a monitoring pro- 
gramme together with H.E.S.S., VERITAS and VLBA (Acciari 
et al. 2009). 

Using the Large Area Telescope (LAT) on-board the Fermi 
Gamma-ray Space Telescope (Fermi-GST) high-energy y- 
radiation above 100 MeV (up to 30 GeV) from M87 was detected 
(Abdo et al. 2009) during the first ten months of the mission. The 
spectral index of 2.26+0.13 is consistent with the VHE measure- 
ments and no significant variability was found. 

The large-scale jet primarily emits from radio to X-rays via 
the synchrotron mechanism. The processes by which the high- 
energy radiation is produced, as well as the regions where these 
processes take place is more debated and has been extensively 
studied in the past years, based especially on radio, optical, 
X-ray, and y-ray observations. The two most promising high- 
energy emission regions are the nucleus and the knot "HST-1". 



At X-rays, Chandra has the angular resolution to separate the 
two components, showing a quite complex behaviour (Harris et 
al. 2009). Unfortunately, the angular resolution of y-ray obser- 
vations does not allow us to resolve the jet structure, but the ob- 
served variability can be used to constrain the size of the emis- 
sion region by requiring that the variability time-scale is longer 
than the light travel time R/c (R is the typical radius of the emis- 
sion region). Even if this information alone is (strictly speaking) 
unable to reveal the location of the emission, it can be corre- 
lated with multiwavelength data in which the source is resolved. 
For the particular case o f M87, studies by Albe rt et al. (2008b)| 
and |Acciari et al. (20 09) have thus favoured the nucleus of M87 
as the most likely origin for the VHE y-radiation (see also 
Abramowski et al. 2012 for a recent summary). 

The discovery of high-energy y-ray emission from M87 (and 
more generally from radio-galaxies) stimulated an intense theo- 
retical work aimed at clarifying the mechanisms that produce the 
observed radiation. As already anticipated, the rapid (~day) vari- 
ability points to a very compact emission region, ruling out early 
models that invoked emission from the large-scale (~kpc) jet 
(Stawarz et al. 2003). Two main theoretical scenarios have been 
discussed in the literature: in the first one the emission is be- 
lieved to arise in the close vicinity of the SMBH through electro- 
magnetic mechanisms acting in the magnetosphere of the black 
hole (e.g. Neronov & Aharonian 2007; Rieger & Aharonian 
2008; Levinson & Rieger 2011). In the other class of models, 
the emission is postulated to originate in the innermost region 
of the jet, as with blazars. The models differ on the geometry 
of the emission regions and on the emission mechanism (e.g. 
Georganopoulos et al. 2005; Lenain et al. 2008; Tavecchio & 
Ghisellini 2008; Giannios et al. 2010; Barkov et al. 2010). In 
the latter class of models, multiple regions are generally as- 
sumed to fully reproduce the observed spectral energy distribu- 
tion (SED). In particular in the structured jet model of Ghisellini 
et al. (2005) and Tavecchio & Ghisellini 2008 one envisages that 
the jet has an inner, faster (bulk Lorentz factor F = 10-15) 
"spine" surrounded by a slower (F = 3-4) sheath. The ra- 
diative interaction between the two regions leads to an effective 
enhancement of the inverse Compton emission and the possibil- 
ity to decelerate the spine through the "Compton drag" effect 
(Begelman & Sikora 1987). This structure may account for the 
emission properties of both blazars (observed at small angles, 
for which the emission is dominated by the spine) and radio- 
galaxies (larger viewing angles and substantial emission from 
the layer). Another possibility is the involvement of electromag- 
netic cascades from relativistic electrons that propagate along 
the jet (Sitarek & Bednarek; and Roustazadeh & Bottcher 2010). 
Finally, other models assume emission from the jet but iden- 
tify the emission region in the HST-1 knot, invoking extreme 
re-collimation of the flow to reproduce the fast variability (e.g. 
Stawarz et al. 2006; Bromberg & Levinson 2009). 

Whereas in the cases mention above the models have been 
applied to interpret high or flaring states of M87, a more detailed 
look into the characteristics of the source's lower emission levels 
is still lacking. In this paper we will focus on the analysis of 
long-term monitoring observations of M87 with MAGIC during 
a low-emission state between 2005-2007. 

2. Observations and data analysis 

The data were collected with the stand-alone MAGIC-I 17 m 
diameter imaging atmospheric Cherenkov telescope situated in 
the Canary Island of La Palma at the Roque de los Muchachos 
Observatory at an altitude of 2200 m above the sea level. The 
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energy threshold of the standard trigger configuration for obser- 
vations at low zenith angles is about 60 Ge V. The angular resolu- 
tion is ~0. 1 ° on an event-to-event basis and the energy resolution 
above 150 GeV is -25% (see Albert et al. 2008a for details). 

MAGIC started to observe M87 regularly in 2005 and con- 
tinued to monitor it during several observing cycles until 2012 
(see e.g. the recent announcement of an enhanced emission state 
in February 2010: Mariotti et al. 2010 and Abramowski et al. 
2012). This paper focuses on data taken between March 2005 
and June 2007 with the stand-alone MAGIC-I telescope. The re- 
sults of the 2008 observing campaign, when M87 entered an in- 
tense high state, have been previously published (see Albert et 
al. 2008b and Acciari et al. 2009 for details). 

The MAGIC-I telescope underwent several significant up- 
grades during the long observation period. These upgrades 
specifically include that optical point spread function (PSF) of 
the 17 m mirror dish was improved due to a better understand- 
ing of the automatic mirror control ( AMC) from about 20 mm 
at the beginning of 2005 to 13 mm at the end of April 2005. In 
the subsequent observing time the telescope underwent several 
distinct mirror-focusing adjustments that resulted in a varying 
optical PSF between 13 mm and 15 mm. While the varying PSF 
has no significant effect on our angular resolution, it does affect 
our collection area, especially near the threshold. These effects 
have been taken into account in the Monte Carlo simulation. 

The readout chain has been significantly enhanced during the 
observing time. While the data recorded until April 2006 used a 
300MSample/s flash analogue-to-digital converter (FADC) sys- 
tem, subsequently a new 2 GSample/s Mux-FADC system was 
installed in February 2007 for the whole MAGIC-I camera. 
During the transition time between April 2006 until February 
2007 a splitter system was installed (and continuously used un- 
til today), which allowed splitting the signal between the Mux- 
FADC readout and the trigger electronics. Accordingly, the data 
discussed in this paper include three different readout stages: 
i) 300MS/s readout without splitters, ii) 300MS/s readout in- 
cluding splitters, and Hi) 2 GS/s readout with splitters. Events 
recorded with the 300MS/s system were stretched up to 6 ns full 
width half maximum (FWHM) (10 ns for the low gain) to en- 
sure proper sampling. With the introduction of the new 2 GS/s 
readout system this stretching as well as the high and low gain 
differentiation became unnecessary. 

MAGIC observations of M87 were performed in two dif- 
ferent modes: ON-OFF, and wobble. In the ON-OFF mode the 
telescope is pointed directly at the source during the ON runs. 
The background is estimated with dedicated OFF runs, which 
are recorded separately in a field of view of the sky where no y- 
ray signal is expected but the night sky background (NSB) and 
zenith angles match the ON observations. It is imperative that 
the telescope conditions are the same for the ON and the OFF 
data to avoid additional systematic effects. Accordingly, only 
OFF data taken during the same observing period as the ON data 
were used in this analysis. Only the 2005 observations of M87 
were recorded in ON-OFF mode (see also Table 1). Data taken 
in the later years were collected in the so-called wobble mode 
(Fomin et al. 1994) where the telescope is not pointed directly 
at the source position, but 0.4° away. Accordingly, wobble data 
do not require dedicated OFF runs since the background can be 
estimated from points in the camera that are equivalent to the 
source position. Those positions are equidistant from the camera 
centre and have sufficient separation from the source location so 
that the expected y-ray signal does not spill over into the OFF 
region. 



Throughout this paper we have applied the same 
timing-based analysis to the entire dataset, as described in 
| Aleksic eT al. (2010a) After applying the standard calibration 
and signal extraction algorithms (Albert et al. 2008c), we 
have performed the so-called "time image cleaning" procedure 
(Aliu et al. 2009), which uses the strictly correlated nature 
of the Cherenkov light pulses (both in time and in space) to 
distinguish from spurious (randomly distributed) signals from 
NSB and electronic noise. The first step requires at least six 
photoelectrons in the core pixels and three photoelectrons in the 
boundary pixels of the images (Fegan 1997). Only pixels with 
at least two neighbouring pixels, whose photons arrive within 
less than 1.75 ns time difference, survive the second cleaning 
step. The third step repeats the cleaning of the second step, but 
requires only one adjacent pixel within the 1.75 ns time window. 

We used the recipe of |Hillas (1985)| to calculate image pa- 
rameters such as WIDTH (the root mean squared of the shower 
image amplitude along the minor axis), LENGTH (the root mean 
squared of the shower image amplitude along the major axis) 
and SIZE (corresponding to the total light content of the im- 
age). The strong domination of air showers induced by charged 
cosmic -rays over y-ray induced showers requires a strong event 
selection. This so-called "y-hadron separation" uses a SIZE- 
dependent parabolic cut in AREA = WIDTH x LENGTH x tt 
(Riegel et al. 2005). The cut parameters were optimised using 
contemporaneous Crab Nebula data. The analysis energy thresh- 
old (defined as the peak of the Monte Carlo energy distribution) 
corresponds to ~ 140 GeV. 

The energy of the primary y-ray events was estimated us- 
ing the Random Forest regression method (Albert et al. 2008d) 
trained on Monte Carlo simulated y-ray events. We used the 
DISP method (Fomin et al. 1994, Lessard et al. 2001) to re- 
construct the arrival direction of the events, which was ex- 
tended to include the shower timing information as described 
in |Aleksic et al. (2010a)l 

Each period with different telescope conditions was first 
analysed separately with matching OFF-data (when needed) and 
Monte Carlo simulations from the corresponding period and fi- 
nally combined after the y-hadron separation. A stacked anal- 
ysis of long-term MAGIC-I observations was first performed 
in |Aleksic et al . (2010b) An excellent agreement between the 
published Crab Nebula spectrum and the stacked analysis was 
found, confirming the stability of the analysis over long time- 
scales. 

The entire observing campaign consists of 154.1 h out of 
which 128.6 h survived the quality selection cuts. This selection 
ensures that only data with good weather conditions and without 
any technical problems are used in the analysis. TableQ]gives an 
overview of how these datasets were distributed over the entire 
campaign. M87 culminates at 16° when observed from La Palma 
and most of the data were taken with low zenith angles (below 
30° with some exceptions up to 45°). 

3. MAGIC results 

The # 2 distribution of the entire observing campaign is shown in 
Fig-El An excess of 862 events over 12454 background events 
is apparent, corresponding to a significance of 7.0 <x (calculated 
with formula 17 of Li & Ma 1983). 

We also investigated the stability of the signal over time. 
Since the observations are not evenly spaced (often interrupted 
by hardware changes of the telescope) and the length of each 
observing run differs greatly, the binning of the light curve is 
not uniform. Instead, we divided the data into calender months, 
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Table 1. List of the effective observing times obtained with each hardware setup. See text for details. 



Obs. period 


Obs. date 


Obs. mode 


read-out system 


PSF 


Obs. time 


2005a 
2005b 


2005-03-15 -2005-04-10 
2005-05-07 - 2005-05-08 


On-Off 
On-Off 


300MS/S 
300MS/S 


20 mm 
13 mm 


6.5 h 
1.5h 


2006a 
2006b 


2005- 12-30 -2006-04-03 

2006- 05-23 -2006-12-31 


Wobble 
Wobble 


300MS/S 
300MS/S + Splitters 


15 mm 
14 mm 


53.4 h 
26.0 h 


2007 


2007-02-13 -2007-05-19 


Wobble 


2 GS/s + Splitters 


14 mm 


42.8 h 



862.2 excess events, 12453.8 background events 
Significance: 7.0 

4500 
4000 
3500 



3000 
2500 
2000 
1500 
1000 
500 



0.05 0.1 0.15 0.2 0.25 0.3 

a 2 [deg 2 ] 

Fig. 1. d 2 distribution of the combined 128.6 h MAGIC observa- 
tions between 2005 and 2007. -& 2 is the squared angular distance 
between the source position and the reconstructed arrival direc- 
tion of each event. The grey shaded area below the red crosses 
is the OFF data-sample, while the black crosses correspond to 
the ON data. The excess is point-like compared to the PSF with 
a significance of 7 standard deviations (calculated with formula 
17ofLi&Ma 1983). 
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Fig. 2. Light curve of the integral y-ray flux above 100 GeV dur- 
ing the 2005-2007 MAGIC observing period. The dashed blue 
line corresponds to the fit result of a linear function to the data 
points, with a reduced x 1 Id.o.f. = 5.66/11 and an average flux 
of (5 .06 + 0.77) x 1 0" 1 2 s" 1 cirT 2 . 



starting from March 2005 until May 2007. Note that these bins 
have been chosen a posteriori and therefore no trial factors are 
involved. We assumed that each positive excess is related to M87 
and derived the corresponding flux value based on the compari- 
son with the Monte Carlo simulation of each period, taking into 
account the collection area, cut efficiency, trigger efficiency, and 
dead-time-corrected observing time. As shown in the light curve 
of Fig. |2] the excess over the entire observing campaign is com- 
patible with a fit to a constant flux of (5.06±0.77)x 10~ 12 s^'cirT 2 
with a reduced x 1 of 0.51 (corresponding to a probability of 
90%). We also analysed each observing night independently, but 
none yielded a significant detection on its own. Fig. [3] shows 
the cumulative excess distribution over time. The excess grows 
linearly, which is consistent with a constant emission, and shows 
that no significant flare occurred during our observing campaign. 

A spectrum was derived using the same analysis technique 
as described for the extraction of the light curve. Additionally, 
a spill-over correction was applied to take into account the fi- 
nite energy resolution (around 25% for the data discussed here). 
Because M87 is a nearby source, the attenuation of VHE y-rays 
by the extra-galactic background light is negligible in our energy 
range (Dommguez et al. 2011). The MAGIC spectrum is well 
described by a single power law: dN/dE = /o (E/300GeV) r 
with a spectral index y = -2.21 + 0.21 and a flux normalization 
f = (7.7 + 1.3)xlO- 8 TeV" 1 s- I m- 2 . The MAGIC spectral points 
and fit are shown in blue in Fig. |4] 

A detailed analysis of the VHE y-ray variability of M87 is 
published elsewhere (Abramowski et al. 2012). The low state 
emission reported in our work is compatible to the flux ob- 
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Fig. 3. Evolution of excess events from M87 over dead-time- 
corrected observing time. The linear increase is consistent with a 
constant flux of the source during the entire observing campaign. 
Only wobble data are shown in this figure in order to reduce sys- 
tematic uncertainties in the ON-OFF subtraction. 



served by the H.E.S.S. Collaboration in 2004 and 2006, as well 
as the 2007 observations with the VERITAS array (see Fig. 1 in 
Abramowski et al. 2012). 

The spectral index derived from our observations (y = 
-2.21 + 0.21) is statistically compatible with previously re- 
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Fig. 4. Combined Fermi- LAT and MAGIC differential energy 
spectrum over four orders of magnitude in energy starting from 
lOOMeV until 2TeV. The Fermi-LAT spectrum was obtained 
from |Abdo et al. (20 09) The blue dashed line corresponds to the 
extrapolation of the fit to the MAGIC data points into the Fermi- 
LAT energy range, while the green line represents the combined 
fit to the MAGIC and Fermi-LAT data. No break or change in 
the spectral slope is apparent and both fits are statistically com- 
patible. Note that the two observations are not contemporaneous, 
but |Abdo et al. (2009)| concluded that the source was most likely 
in a low-emission state during the observations. 



ported results from H.E.S.S. (y=2.62+0.35 for the 2004 data 
and y=2.22+0.15 for 2005, Aharonian et al. 2006) as well as 
VERITAS (y=2.31+0.17 in 2007, Acciari et al. 2008). 

During the MAGIC observing periods no observations were 
performed in the GeV energy band. However, for comparison, 
we report (red squares, Fig. 4) the Fermi-LAT spectrum ob- 
tained by integrating over the first ten months of all-sky survey 
data (Aug. 4, 2008 - May 31, 2009, Abdo et al. 2009). Based on 
contemporaneous X-ray and radio data, |Abdo et al. (2009)] ar- 
gued that M87 was most likely in a low-emission state during 
the Fermi-LAT observations. To investigate this hypothesis in 
more detail, we performed a combined fit to the MAGIC and 
Fermi-LAT results (green line, Fig. 4). The combined fit yields 
a spectral index y - -2.17 ± 0.03 and a flux normalization at 
300GeV f = (7.1 + 1.0) x 10~ 8 Te^'s^nT 2 . The reduced x 2 
changes from 1.14 to 0.86. The fit result is statistically compati- 
ble with the independent fit to the MAGIC data and confirms the 
low state during the Fermi-LAT observations. There is no indi- 
cation of a break or change in the spectral slope. We therefore 
included the Fermi-LAT data in our modelling of the spectral 
energy distribution in the next section. It is worth noting that the 
peak of the spectral energy distribution appears to lie at particu- 
larly low energies, below 100 Me V. 

A detailed discussion of the potential systematics can be 
found in Alber fet al. (2008a)| As a summary, the systematic en- 
ergy scale error is about 16%, the error of the flux normalization 
is estimated to be 11% and the systematic slope error is ±0.2. 
Additionally, we have to consider how the combination of dif- 
ferent observing modes and telescope configurations affects the 
systematic errors. We have combined Monte Carlo simulations 
that were tailored for each observing epoch (including the re- 
spective PSF) and tried to reconstruct the original Monte Carlo 
spectrum, which was fixed to a power law with a spectral index 



of -2.2 (identical to the detected spectral index of M87). From 
the maximum error in the reconstructed Monte Carlo spectrum 
we find that the error of the flux normalization should be <10% 
and the error of the slope should be <0.03 (note that this applies 
only to hard spectra with a spectral index <2.5, as is the case 
for the data discussed here). Finally, an independent second data 
analysis yielded compatible results. 

4. Modelling the SED 

Fig. |5] shows the SED of the core of M87, adapted from 
Tavecchio & Ghisellini (2008) Green open symbols represent 
historical data, black filled circles show the MAGIC spectrum. 
For comparison we also show the H.E.S.S. spectra taken in 2004 
and 2005 (magenta triangles and red open circles, respectively, 
from Aharonian et al. 2006) and the 2007 VERITAS spectrum 
(cyan squares, from Acciari et al. 2008). The 2005-2007 spec- 
trum was very similar to that measured by H.E.S.S. in the "low" 
2004 state and by VERITAS in 2007. The green "bow-tie" re- 
ports the nuclear X-ray spectrum measured by Chandra in 2000 
(Balmaverde et al. 2006). The black "bow-tie" shows the average 
X-ray spectrum during the period covered by the MAGIC obser- 
vations, obtained by scaling the spectrum from the year 2000 
by a factor of four, as inferred from the light curve reported 
by |Harris et al. (2009)| This procedure was necessary since no 
spectrum is available for the period of the MAGIC observations 
and both campaigns were not organised jointly, which results 
in very little overlapping observing time. The stability of the 
low-state spectrum supports this procedure of scaling. As ex- 
plained in Section 3, we also included the non-contemporaneous 
Fermi-LAT spectrum (orange open diamonds) obtained by inte- 
grating over the first ten months of all-sky survey data (Abdo et 
al. 2009). The SED shows two pronounced bumps, one peaking 
in the IR band, one extending from MeV to TeV energies. 

As discussed in the introduction, models of the VHE emis- 
sion from M87 include a class that considers the emission from 
ultra-relativistic electrons accelerated in the black hole magne- 
tosphere in the vicinity of the horizon, and another group as- 
suming emission from the jet. In the latter group some scenar- 
ios are based on variations of the general framework adopted 
for blazars, i.e. assuming that the emission comes from the 
innermost regions of the jet, corresponding to the radio core 
(at 100-1000 Schwarzschild radii from the SMBH), while 
others propose that the radiation is produced at larger dis- 
tances, possibly at a shock generated by a strong re-collimation 
of the flow. In the context of jet models, as discussed in 
Tavecchio & Ghisellini (2008), one-zone models face a severe 
problem in describing the whole SED of M87, from radio up 
to TeV energies. Indeed, the nuclear environment of M87 does 
not show evidence of significant sources of soft photons (e.g. 
from the torus or star formation sites) and therefore the in- 
verse Compton emission is likely dominated by the SSC com- 
ponent. In that case, it is possible to derive the value of the 
required Doppler factor, 6, which depends only on the syn- 
chrotron and SSC peak frequencies and luminosities and on 
the variability timescale. For an SED such as the one dis- 
played by M87, whose synchrotron and SSC components peak 
in the IR and TeV band, respectively, one infers unacceptably 
high values of the Doppler factor, 6 > 100. The most di- 
rect solution to this problem is to assume that the jet emis- 
sion region is structured, as in the "decelerating jet" model 
of Georganopoulos et al. (2005) or in the "spine-layer" model 
of Ghisellini, Tavecchio & Chiaberge (2005) Other possibilities 



include multiple regions moving into a wider jet, as expected 
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from the inner regions before the initial collimation (Lenain et 
al. 2008), or emission by several randomly oriented active re- 
gions resulting from reconnection events in the jet (the so-called 
"jets-in-the-jet" scenario of Giannios et al. 2009; Giannios et al. 
2010). By construction, all these models predict that M87, if ob- 
served at smaller angles, would display an SED resembling those 
commonly observed for blazars. 

We modelled the SED using the structured-jet model of 
Ghisellini, Tavecchio & Chiaberge (2005) previously applied to 
M87 in |Tavecchio & Ghisellini (2008)| assuming that the low- 
state y-ray emission is coming from the radio core, as is sus- 
pected for the high-state emission (Acciari et al. 2009). The 
model assumes that the jet has an inner fast core (spine) with 
a bulk Lorentz factor F$, surrounded by a slower layer, with 
a bulk Lorentz factor Fl- In both regions, relativistic elec- 
trons emit through synchrotron and inverse Compton mecha- 
nisms. The velocity structure plays an important role in deter- 
mining the emission properties of the jet. Indeed, the radia- 
tive interaction between the layer and the spine results in the 
amplification of the inverse Compton emission of both com- 
ponents. In the rest frame of each component the emission of 
the other is amplified because of the relative speed within the 
two regions: this "external" radiation contributes to the total 
energy density, enhancing the emitted inverse Compton radi- 
ation. Depending on the parameters, this "external Compton" 
emission can dominate the internal synchrotron self-Compton 
(SSC) component that, especially in TeV blazars, is suppressed 
because scatterings mainly occur in the Klein-Nishina (KN) 
regime. We refer to Ghisellini, Tavecchio & Chiaberge (2005) 
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and |Tavecchio & Ghisellini (2008)| for a full description of the 
model. 

The model is completely specified by the following parame- 
ters: i) the spine is assumed to be a cylinder of radius R, height 

(as measured in the spine frame) and in motion with bulk 
Lorentz factor Fs; ii) the layer is modelled as a hollow cylinder 
with internal radius R, external radius Ri, height (as mea- 
sured in the frame of the layer), and bulk Lorentz factor Fl- Each 
region contains a tangled magnetic field with intensity Z?s and 
Z?l respectively, and is filled by relativistic electrons assumed 
to follow a (purely phenomenological) smoothed broken power- 
law distribution extending from y m i„ to y max with indices n\, 
below and above the break at yt,. The normalization of this distri- 
bution is calculated assuming that the system produces a (bolo- 
metric) synchrotron luminosity L syn (as measured in the local 
frame), which is an input parameter of the model. We assume 
that //l > Hs. As noted above, the seed photons for the IC scat- 
tering are not only those produced locally in the spine (layer), but 
we also consider the photons produced in the layer (spine). The 
result of the modelling is reported in Fig. [5] where we show the 
SED of the emission produced by the spine (red) and the layer 
(blue) and their sum (black). The adopted parameters are re- 
ported in Table [2] The parameters we adopted are very similar to 
those derived in Tavecchio & Ghisellini (2008) for the high 2005 
state. In particular, we used the same values of the spine and 
layer dimensions and Doppler factors. Our model is therefore 
compatible with our initial assumption that the low- and high- 
state emissions are originating from the same emission region 
near the radio core. Since in the SED considered in Tavecchio 
& Ghisellini (2008) there were no measures in the GeV band, 
the model of the flux in that band was unconstrained. The LAT 
data we used in the present model instead allow us to constrain 
the inverse Compton bump of the spine to a luminosity slightly 
lower than that assumed in Tavecchio & Ghisellini (2008): this 
difference accounts for a stronger magnetic field derived here for 
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Fig. 5. SED of the core of M87 (green open squares) together 
with the MAGIC data points (black filled circles). For compari- 
son we report the H.E.S.S. spectra taken in 2004 (magenta tan- 
gles) and 2005 (red open circles), from |Aharonian et~a l. (2006) 
and the VERITAS spectrum taken in 2007 (Acciari et al. 2008, 
cyan filled squares). The green bow-tie reports the X-ray spec- 
trum as measured by Chandra in 2000 (from Balmaverde et al. 
2006). As in |Tavecchio & Ghisellini (2008)| we reproduced the 
X-ray emission of the period 2005-2007 assuming the same 
slope and a higher normalization (black bow-tie). In the high- 
energy y-ray range the green open square corresponds to the 
EGRET upper limit and the orange diamonds to the Fermi-LAH 
energy spectrum from Abdo et al. 2009. The lines report the 
emission from the spine (red) and from the layer (blue), cal- 
culated with the parameters reported in Table 2, and their sum 
(black). 



the spine. Similarly to Tavecchio & Ghisellini (2008), since we 
are forced to reproduce a high X-ray state, the model inevitably 
over-predicts the non-simultaneous flux in the IR-optical bands. 
However, as already noted in Tavecchio & Ghisellini (2008), be- 
cause the X-ray and optical fluxes are correlated (Perlman et al. 
2003), we expect that high X-ray fluxes also corresponds to op- 
tical states higher than those presented in the SED. 

A feature of the "spine-layer" model for M87 is that the 
IC emission from the layer, which accounts for the observed 
VHE emission, is partly absorbed through the interaction with 
the optical-IR photons produced in the spine and the subsequent 
production of electron positron pairs. This leads to a relatively 
soft spectrum, suitable for reproducing the low-level spectrum 
measured by MAGIC (which extends to ~2TeV), but difficult 
to reconcile with the hard spectrum recorded during the high 
state by H.E.S.S. A possible solution of this problem is to en- 
large the emission regions, since this would reduce the den- 
sity of the target photons. However, the increase of the source 
size is limited by the observed short variability time-scales, 
R < 5 x IQ 15 6 cm. In summary, we are able to describe the VHE 
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Table 2. Input parameters of the models for the layer and the spine shown in Fig. [5] All quantities (except the bulk Lorentz factors F 
and the viewing angle 0) are measured in the rest frame of the emitting plasma. The external radius of the layer is fixed to the value 
R2 = 1.2 x R. A detailed description of the model and its parameters can be found in the text. 
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low-state emission with parameters similar to those adopted 
in |Tavecchio & Ghisellini (2008)| to reproduce the high state. 
We can therefore conclude that in the framework of the model 
adopted here, both high and quiescent states can originate from 
the same emitting region. The opacity of the emitting region to 
the y-ray photons is a difficulty that also affects other models, 
especially those focusing on the acceleration in the BH magne- 
tosphere, where the environment is expected to be extremely rich 
in IR photons that originate in the accretion flow (e.g. Neronov & 
Aharonian 2007; Li et al. 2009). This problem, instead, could be 
relaxed in the "mini-jets" model of Giann ios et al. (2010)| due 
to the high Doppler factor 5^10 assumed to characterize the 
emitting regions. We also note that in the structured-jet scenario 
presented here, a strict correlation between these two spectral 
components is not required (even if it is possible), because the 
MeV-GeV emission originates in the spine while the VHE radi- 
ation is produced in the layer. 

Another point concerns the jet power required to reproduce 
the SED. The inferred jet power (dominated by the Poynting 
flux associated to the magnetic field) is Pj et = 1.5 x 10 44 erg 
s _1 . The accretion luminosity in M87 is estimated to be about 
L accr = 10 40 - 10 41 erg s _1 (e.g. Perlman et al. 2001), which for 
a black hole mass of M — 6 x 10 9 M corresponds to Eddington 
ratios of L a cci7-^Edd = 10 8 - 10~ 7 . The luminosities as low as 
this, common among low -power radio galaxies are currently in- 
terpreted to be the result of a radiatively inefficient flow, with ef- 
ficiencies of about rj = L aca /Mc 2 = 10 4 — 10 3 (e.g. Balmaverde 
et al. 2008). In that case the power request of the jet can be easily 
fulfilled by the accretion power, Mc 2 = L aca /j] ~ 10 45 erg s" 1 . 
Finally, we also note that a magnetically dominated jet (Poynting 
to kinetic power ratio higher than ~100) is postulated to power 
the "mini-jets" in the |Giannios et al. (2010)| model. 

A final comment concerns the angle within the plasma veloc- 
ity and the line of sight, which in our model we assumed to be 
0=18 deg. This value is well within the range confidently de- 
rived through radio observations, 6 - 15-25 deg (Acciari et al. 
2009). Angles larger then ~ 20 deg would imply a de-boosting of 
the spine emission (see Fig. 2 of Tavecchio & Ghisellini 2008), 
with the consequent increase in the intrinsic luminosity of the 
spine. The increased luminosity, in turn, would imply a higher 
power of the jet and, more importantly, would increase the opti- 
cal depth to VHE photons emitted by the layer. 

5. Summary 

MAGIC has detected weak and steady VHE y-ray emission from 
M87 between 2005 and 2007. The flux and spectral shape are 
consistent with a straight power law extrapolation of the pub- 
lished Fermi-LAT spectrum at lower energies although both ob- 
servations are not contemporaneous. Our measurements are also 
compatible with previously reported low states by VERITAS and 
H.E.S.S., which in turn suggests that the observed low-emission 
level and spectral characteristics are stable over a long time pe- 
riod. 



We were able to describe this emission with a structured jet 
model, which separates the jet into a spine and an outer layer. 
The assumption that the low-state VHE y-ray emission origi- 
nates from the radio core of the jet like the high-state emission 
is consistent with our model. It should be noted, however, that 
the low-state emission alone cannot constrain the emission re- 
gion due to the lack of variability of the measured signal and the 
softer VHE y-ray spectrum. The parameters we derived from 
the model fit are consistent (even identical for the spine/layer 
dimensions and Doppler factors) with Tavecchio & Ghisellini 
(2008) for the high 2005 state. Because the present model in- 
cludes GeV data from the Fermi-LAT, we were able to constrain 
the inverse Compton bump of the spine to a slightly lower lumi- 
nosity compared to the assumptions in Tavecchio & Ghisellini 
(2008), which results in a stronger magnetic field for the spine. 

An interesting feature of the model is the transition region 
between the emission from the spine and that of the layer be- 
tween 40 and 100 GeV (the exact position depends on the cho- 
sen model parameters, see Fig. [5] and Table 2). Currently, no 
measurements of M87 are available in this energy range and the 
non-contemporaneity of the multiwavelength data limits the ac- 
curacy of the spine-layer model. Since fall 2009, MAGIC con- 
sists of two 17 m diameter telescopes, which observe in a stereo- 
scopic mode. This upgrade has significantly improved the in- 
strument's capabilities at energies below 100 GeV as reported in 
Aleksic et al. (2012), It is therefore possible that a future, deep 
observation of M87 with the stereoscopic MAGIC system and 
contemporaneous multiwavelength data will reveal a feature in 
the otherwise smooth power law spectrum of M87, confirming 
or disproving the validity of the structured jet model. 
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